Leak-rate of seals: comparison of theory with experiment 
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Seals are extremely useful devices to prevent fluid leakage. We present experimental results for the 
leak-rate of rubber seals, and compare the results to a novel theory, which is based on percolation 
theory and a recently developed contact mechanics theory. We find good agreement between theory 
and experiment. 



A seal is a device for closing a gap or making a joint 
fluid tight [ij. Seals play a crucial role in many modern 
engineering devices, and the failure of seals may result in 
catastrophic events, such as the Challenger disaster. In 
spite of its apparent simplicity, it is not easy to predict 
the leak-rate and (for dynamic seals) the friction forces [2] 
for seals. The main problem is the influence of surface 
roughness on the contact mechanics at the seal-substrate 
interface. Most surfaces of engineering interest have sur- 
face roughness on a wide range of length scales [3], e.g, 
from cm to nm, which will influence the leak rate and 
friction of seals, and accounting for the whole range of 
surface roughness is impossible using standard numerical 
methods, such as the Finite Element Method. 

In this paper we present experimental results for the 
leak-rate of rubber seals, and compare the results to 
a novel theory [1, 0, [Hf, which is based on percola- 
tion theory and a recently developed contact mechan- 
ics theory la i, s, 113, m m, which accurately takes 

into account the elastic coupling between the contact re- 
gions in the nominal rubber-substrate contact area. Ear- 
lier contact mechanics models, such as the Greenwood- 
Williamson [13] model or the model of Bush et al[3], ne- 
glect this elastic coupling, which results in highly incor- 
rect results [15|,[l6|, in particular for the relations between 
the squeezing pressure and the interfacial separation[l7|. 
We assume that purely elastic deformation occurs in the 
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FIG. 1: Rubber seal (schematic). The liquid on the left-hand- 
side is under the hydrostatic pressure Pa and the liquid to the 
right under the pressure Pb (usually, Pb is the atmospheric 
pressure). The pressure difference AP = Pa — Pb results in 
liquid flow at the interface between the rubber seal and the 
rough substrate surface. The volume of liquid flow per unit 
time is denoted by Q, and depends on the squeezing pressure 
Po acting on the rubber seal. 
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FIG. 2: The contact region at different magnifications C = 3, 
9, 12 and 648, are shown in (a)-(d) respectively. When the 
magnification increases from 9 to 12 the non-contact region 
percolate. At the lowest magnification C = l' A(l) = Aq. The 
figure is the result of Molecular Dynamics simulations of the 
contact between elastic solids with randomly rough surfaces, 
see Ref. [5]. 



solids, which is the case for rubber seals. 

Consider the fluid leakage through a rubber seal, from 
a high fluid pressure Pa region, to a low fluid pressure Pb 
region, as in Fig. [TJ Assume that the nominal contact 
region between the rubber and the hard countersurface 
is rectangular with area x Ly. We assume that the 
high pressure fluid region is for x < and the low pres- 
sure region for x > L^. We "divide" the contact region 
into squares with the side Lx = L and the area Aq = 
(this assumes that N = Ly/L^ is an integer, but this 
restriction does not affect the final result). Now, let us 
study the contact between the two solids within one of 
the squares as we change the magnification C,. We define 
C = L/\^ where A is the resolution. We study how the 



apparent contact area (projected on the x?/-plane), 
between the two sohds depends on the magnification 
At the lowest magnification we cannot observe any sur- 
face roughness, and the contact between the sohds ap- 
pears to be complete i.e., A{1) = Aq. As we increase the 
magnification we will observe some inter facial roughness, 
and the (apparent) contact area will decrease. At high 
enough magnification, say C = Cc a percolating path of 
non-contact area will be observed for the first time, see 
Fig. [21 We denote the most narrow (and least high) con- 
striction along this percolation path as the critical con- 
striction. The critical constriction will have the lateral 
size Ac = L/(c and the surface separation at this point 
is denoted by Uc. We can calculate Uc ~ ui{CX using a 
recently developed contact mechanics theory Thus, 
we define ui{() to be the (average) height separating the 
surfaces which appear to come into contact when the 
magnification decreases from ( to ( — where A( is a 
small (infinitesimal) change in the magnification. ui{Q 
can be calculated as described below. As we continue to 
increase the magnification we will find more percolating 
channels between the surfaces, but these will have more 
narrow constrictions than the first channel which appears 
at C = Cc, and as a first approximation we will neglect 
the contribution to the leak-rate from these channels [5]. 

A first rough estimate of the leak-rate is obtained by 
assuming that all the leakage occurs through the critical 
percolation channel, and that the whole pressure drop 
AP = Pa — Pb (where Pa and Pb is the pressure to the left 
and right of the seal) occurs over the critical constriction 
[of width and length Xc ^ L/(c and height Uc = ui{(c)]' 
If we approximate the critical constriction as a pore with 
rectangular cross section (width and length Ac and height 
Uc << Ac), and if assume an incompressible Newtonian 
fiuid, the volume-fiow per unit time through the critical 
constriction will be given by (Poiseuille fiow) 



perpendicular to the fiuid fiow. In addition, the chan- 
nel is of course not exactly rectangular in the x?/-plane, 
and this too will effect a. Note also that a given per- 
colation channel could have several narrow (critical or 
nearly critical) constrictions of nearly the same dimen- 
sion which would reduce the fiow along the channel. But 
in this case one would also expect more channels from the 
high to the low fiuid pressure side of the junction, which 
would tend to increase the leak rate. These two effects 
will, at least in the simplest picture, compensate each 
other (see discussion in Ref. [5]). Finally, since there are 
N = Ly/Lx square areas in the rubber-countersurface 
(apparent) contact area, we get the total leak-rate 
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To complete the theory we must calculate the separa- 
tion Uc = Ui{(c) of the surfaces at the critical constric- 
tion. We first determine the critical magnification by 
assuming that the apparent relative contact area at this 
point is given by site percolation theory. Thus, the rel- 
ative contact area A{Q/Ao ^ 1 — where Pc is the so 
called site percolation threshold[18]. For an infinite-sized 
systems Pc ~ 0.696 for a hexagonal lattice and 0.593 for 
a square lattice[18]. For finite sized systems the perco- 
lation will, on the average, occur for (slightly) smaller 
values of p, and fiuctuations in the percolation threshold 
will occur between different realizations of the same phys- 
ical system. We take Pc ^ 0.6 so that A{(c)/Ao ~ 0.4 
will determine the critical magnification C = Cc- 

The (apparent) relative contact area A{Q/Ao at the 
magnification ( can be obtained using the contact me- 
chanics formalism developed elsewhere |6[ [sl, [qI, [io|, El, 
where the system is studied at different magnifications 
We havelSl, !?! 
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where r] is the fiuid viscosity. In deriving (1) we have 
assumed laminar fiow and that Uc << Ac, which is al- 
ways satisfied in practice. We have also assumed no-slip 
boundary condition on the solid walls. This assumption 
is not always satisfied at the micro or nano-scale, but is 
likely to be a very good approximation in the present case 
owing to surface roughness which occurs at length-scales 
shorter than the size of the critical constriction. 

In (1) we have introduced a factor a which depends 
on the exact shape of the critical constriction, but which 
is expected to be of order unity. The fiow rate expected 
for a channel with rectangular cross section (height ui 
and width and length Ac with ui « Ac) correspond 
to a = 1. However, the actual fiow channel will not 
have a rectangular cross section but the pore height must 
go continuously to zero at the "edges" in the direction 
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where the surface roughness power spectrum 
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where (...) stands for ensemble average. Here E and u 
are the Young's elastic modulus and the Poisson ratio of 
the rubber. The height profile h{x.) of the rough surface 
can be measured routinely today on all relevant length 
scales using optical and stylus experiments. 

We define ui{Q to be the (average) height separating 
the surfaces which appear to come into contact when the 
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FIG. 3: An asperity contact region observed at the magnifica- 
tion It appears that complete contact occur in the asperity 
contact region, but when the magnification is increasing to the 
highest (atomic scale) magnification (i , it is observed that the 
solids are actually separated by the average distance u{(). 



magnification decreases from ( to ( — where is a 
small (infinitesimal) change in the magnification. ui{() 
is a monotonically decreasing function of and can be 
calculated from the average interfacial separation u{Q 
and A{() using (see Ref. [11]) 

ui{o = u{o+u'm{o/A'{o- 

The quantity u{Q is the average separation between the 
surfaces in the apparent contact regions observed at the 
magnification see Fig.O It can be calculated from|ll| 
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where p{C) = Po^oM(C) and 



nq \ -1/2 

tt/ dq'q'^dq') 



The function P{q,p, C) is given by 
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where s{qX) =w{qX)/E\ 

We have performed a very simple experiment to test 
the theory presented above. In Fig. [H we show our set- 
up for measuring the leak-rate of seals. A glass (or 
PMMA) cylinder with a rubber ring (with rectangular 
cross-section) attached to one end is squeezed against a 
hard substrate with well-defined surface roughness. The 
cylinder is filled with water, and the leak-rate of the fluid 
at the rubber-countersurface is detected by the change 
in the height of the fluid in the cylinder. In this case 
the pressure difference AP = — = pgH, where 
g is the gravitation constant, p the fluid density and 
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FIG. 4: Experimental set-up for measuring the leak-rate 
of seals. A glass (or PMMA) cylinder with a rubber ring 
attached to one end is squeezed against a hard substrate 
with well-defined surface roughness. The cylinder is filled 
with water, and the leak-rate of the water at the rubber- 
countersurface is detected by the change in the height of the 
water in the cylinder. 
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FIG. 5: Surface roughness power spectrum of sand paper 120. 
The surface has the root- mean- square roughness 44 /xm and 
the surface area (including only the surface roughness with 
wavelength above Ai = 20 /xm) is about 40% larger than the 
nominal surface area Ao (i.e., the surface area projected on 
the x^-plane). 



H the height of the fluid column. With ^ 1 m 
we get typically AP ^ 0.01 MPa. With the diameter 
of the glass cylinder of order a few cm, the condition 
Po >> AP (which is necessary in order to be able to 
neglect the influence on the contact mechanics from the 
fluid pressure at the rubber-countersurface) is satisfied 
already for loads (at the upper surface of the cylinder) 
of order kg. In our study we use a rubber ring with 
the Young's elastic modulus E = 2.3 MPa, and with 
the inner and outer diameter 3 cm and 4 cm, respec- 
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FIG. 6: Square symbols: the measured leak rate for different 
fluid pressure drop AP = — Pb for the nominal squeezing 
pressure Po ~ 60 kPa. 

tively, and the height 0.5 cm. The rubber ring was 
made from a silicon elastomer (PDMS) prepared using 
a two-component kit (Sylgard 184) purchased from Dow 
Corning (Midland, MI). The kit consist of a base (vinyl- 
terminated polydimethylsiloxane) and a curing agent 
(methylhydrosiloxane-dimethylsiloxane copolymer) with 
a suitable catalyst. From these two components we pre- 
pared a mixture 10:1 (base/cross linker) in weight. The 
mixture was degassed to remove the trapped air induced 
by stirring from the mixing process and then poured into 
casts. The bottom of these casts was made from glass to 
obtain smooth surfaces. The samples were cured in an 
oven at 80° C for 12 h. The substrate is a corundum pa- 
per (grit size 120) with the root-mean-square roughness 
44 /im. From the measured surface topography we ob- 
tain the surface roughness power spectrum C{q) shown 
in Fig. El 

According to (1) we expect the leak-rate to increase 
linearly with the fluid pressure difference AP = Pa — Pb- 
We first performed some experiments to test this pre- 
diction. In Fig. [6] we show the measured leak rate for 
different fluid pressure drop AP for the nominal squeez- 
ing pressure Pq ~ 60 kPa. To within the accuracy of the 
experiment, the leak-rate depends linearly on AP. 

In Fig. Owe show the measured leak rate for ten dif- 
ferent squeezing pressures (square symbols). The solid 
line is the calculated leak rate using the measured rub- 
ber elastic modulus £^ = 2.3 MPa and the surface power 
spectrum C{q) shown in Fig. [H In Fig. [8] we show the 
same as in Fig. [3 but for a logarithmic (with 10 as basis) 
leak-rate scale. 

In Fig. [9] we show the calculated critical pore size as a 
function of the squeezing pressure. Note that the height 
of the critical pore is about 10 times smaller than the 
lateral size of the pore. Finally, in Fig. [10] we show 
the critical magnification where the non-contact area 
percolate, as a function of the squeezing pressure. Note 



FIG. 7: Square symbols: the measured leak rate for ten dif- 
ferent squeezing pressures. The experiment was performed 
twice, corresponding to the two data points for each pres- 
sure. Solid line: the calculated leak rate using the mea- 
sured surface topography, the measured rubber elastic mod- 
ulus P = 2.3 MPa and the fluid pressure difference AP = 
Pa — Pb = 10 kPa obtained from the height of the water col- 
umn. We have used a = 0.2. 
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FIG. 8: The same as in Fig. [71 but for a logarithmic (with 10 
as basis) leak-rate scale. 

that, as expected, the percolation of the non-contact area 
occur at higher and higher magnification as the squeezing 
pressure increases. 

Sand paper has much sharper and larger roughness 
than the counter surface used in normal rubber seal ap- 
plications. However, from a theory point of view it should 
not really matter on which length scale the roughness oc- 
curs, except for "complications" such as the influence of 
adhesion and fluid contamination particles (which tend 
to clog the flow channels). Nevertheless, the theory as- 
sumes that the average surface slope is not too large and 
we plan to measure the leak rate for rubber seal in con- 
tact with sand blasted Plexiglas with a root-mean-square 
roughness in the micrometer range. Our initial experi- 
ment with Plexiglas showed that the leak rate decreased 
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FIG. 9: Calculated critical pore size as a function of the 
squeezing pressure. For the same system as in Fig. [T] 
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FIG. 10: Calculated critical magnification ("c, where the non- 
contact area percolate, as a function of the squeezing pressure. 
For the same system as in Fig. [T] 

by time and finally no leaking could be observed. But 
this experiment used unfiltered tap water which contains 
contamination particles which clogged the channels. We 
are now using distilled water and find the leak rate (for 
a given fluid pressure difference) to be practically time 
independent. But these are still preliminary studies and 
we will report on the final results elsewhere. 

Finally, we note that it is nearly impossible to calcu- 
late a theoretically and we have fitted a to reproduce 
the experimental data. However, this is just one param- 
eter and the dependence of the leak rate on the nomi- 
nal rubber-countersurface pressure is highly non-trivial 
and accurately given by the theory. Also, a is of order 
unity as expected from theory. The only way (as we see 
it) to obtain a theory-estimate of a would be to gener- 
ate critical constrictions using numerically exact contact 
mechanics calculations (e.g., molecular dynamics), and 
to simulate the fluid flow through the constrictions using 
numerical methods of fluid flow dynamics. 



To summarize, we have compared experimental data 
with theory for the leak-rate of seals. The theory is based 
on percolation theory and a recently developed contact 
mechanics theory. The experiments are for silicon rub- 
ber seals in contact with sand paper. The elastic prop- 
erties of the rubber and the surface topography of the 
sand paper are fully characterized. The calculated leak- 
rate Q is in good agreement with experiment. The the- 
ory only account for fluid flow through the percolation 
channels observed at (or close to) the percolation thresh- 
old. A more accurate treatment should include also flow 
channels observed at higher magnification. This prob- 
lem has similarities to current flow in random resistor 
networks pj, 
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